Electromagnetic energy within a magnetic infinite cylinder and scattering properties 

for oblique incidence 
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In this work we analytically calculate the time-averaged electromagnetic energy stored inside 
a nondispersive magnetic isotropic cylinder which is obliquely irradiated by an electromagnetic 
plane wave. An expression for the optical-absorption efficiency in terms of the magnetic internal 
coefficients is also obtained. In the low absorption limit, we derive a relation between the normalized 
internal energy and the optical-absorption efficiency which is not affected by the magnetism and 
the incidence angle. This mentioned relation, indeed, seems to be independent on the shape of 
the scatterer. This universal aspect of the internal energy is connected to the transport velocity 
and consequently to the diffusion coefficient in the multiple scattering regime. Magnetism favors 
high internal energy for low size parameter cylinders, which leads to a low diffusion coefficient for 
electromagnetic propagation in 2D random media. 

PACS numbers: 



I. INTRODUCTION 

The interest in the study of magnetic materials and 
their optical properties has been increased in the ap- 
plied science. Recently, applications on magnetic 2D and 
3D photonic band gaps (PBGs) [JQ, microwave filters, 
metamaterials @ , high density magnetic recording media 
@ , weak localization of light HI have been reported. 
All these works explore the fact that, at microwave or 
radio frequencies, the magnetic materials exhibit large 
values of magnetic permeability 0, 0] ■ 

Electromagnetic (EM) scatteringby magnetic spheres 
has been applied by Kerker et al. [9j. Some unusual fea- 
tures present in single and multiple magnetic Mie scat- 
tering, such as forward-backward asymmetry with prefer- 
ential backward scattering and resonance effects [3, El, 
and vanishing of the energy-transport velocity even for 
small size parameters 0,11], have been studied. In a re- 
cent paper, we have calculated the EM energy stored in- 
side a magnetic sphere and we have shown that, even for 
size parameters much smaller than unity (Rayleigh size 
region), this quantity is strongly enhanced, with sharp 
resonance peaks [l2| . 

The problem of EM scattering by an isotropic circu- 
lar cylinder is not new [l3|, Il4j . For a general case of 
oblique incidence and magnetic scatterers, an analytical 
solution has been provided a long time ago by Wait 15 1 
and has also been treated by Lind and Greenberg [16 1 
in the context of dielectric infinite cylinders. Although 
expressions for the stored energy in a normally illumi- 
nated dispersive and nondispersive dielectric cylinders 
have been documented in the literature [l7| . no analyt- 
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ical study of the magnetism influence at a general case 
of oblique incidence has been so far performed. Our aim 
is to fill this gap with a detailed study of a cylindrical 
magnetic scatterer illuminated at an arbitrary incidence 
angle. We devote special attention to the fields inside the 
scattering center and their application on the calculation 
of the energy-transport velocity in a disordered magnetic 
medium [ll []J|. 

The framework in which the scattering quantities are 
calculated is presented in Sec. |TT]of this description. Es- 
sentially, we present the main expressions obtained solv- 
ing the macroscopic Maxwell's equations for the EM in- 
ternal fields [H, EH . We have adopted the same notation 
as Bohren and Huffman [Hj. In Sec. EH for a general 
case of oblique incidence, we use new relations among 
Bessel function to calculate the scattering coefficients and 
consequently the normalized average EM energy stored 
inside a magnetic infinitely long cylinder. This extends 
the study of [13, H3| for nondispersive scatterers. Also, 
in Sec. IIV1 we derive exact and approximated expres- 
sions for the optical-absorption efficiency in terms of the 
internal coefficients. This last result, which has an ana- 
logue in the single Mie scattering [12j , is important to link 
measurable quantities with the time-averaged EM energy 
[HI, [20|. Numerical results are shown in Sec.lVl Specially, 
for a weakly absorptive magnetic cylinder, we determine 
a relation between the internal energy-enhancement fac- 
tor and the absorption efficiency, which does not depend 
on the incidence angle and the polarization of the incident 
EM wave. Comparing our result with that one obtained 
firstly by [20|, we ascribe the achieved difference only to 
the cylindrical and spherical geometries. This geomet- 
rical consideration allows us to write a relation which 
is i ndep endent on the shape of the scatterer. Following 
[l7l fl8j . we present an application of our calculations to 
the study of the energy-transport velocity. For a two- 
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dimensional disordered magnetic medium, we obtain an 
oscillatory behavior of the energy-transport velocity, as 
a function of the size parameter, even in the Rayleigh 
size region. Briefly, in Appendix [XJ we present approxi- 
mations for the far-field scattering coefficients and some 
considerations about the degree of polarization for mag- 
netic cylinders at normal incidence. 



II. BASIC THEORY 

Let the scatterer be an infinite right circular cylin- 
der with finite radius a embedded in an infinite non- 
absorptive medium. Both the cylinder and the sur- 
rounding medium are assumed to be linear, homoge- 
neous and isotropic, with inductive capacities (ei, /ii) and 
(e, /i), respectively. The incident EM wave, which inter- 
acts with the cylinder, is a plane and monochromatic 
complex wave, with time-harmonic dependence given by 
exp(— iu)t) [HI, [2l|. The quantity w is the angular fre- 
quency and it is considered to be the same for the incident 
and scattered waves (elastic scattering). In addition, sup- 
pose that these media are electromagnetically source- free 
and adopt the international system of units. On account 
of symmetry, the cylindrical scatterer imposes two basic 
linear polarizations for the incident EM wave [III [22j]. 
They are referred to as the TM (or case I) and TE (or 
case II) modes [H|. In the former, the incident electric 
field is parallel to the xz plane, while in the latter it is 
perpendicular to this plane (l9l [2l| - |23j . For both cases, 
consider that £ is the angle between the wavevector k and 
the z axis, where k = |k| = u^/xe) 1 / 2 is the wavenumber. 

Inside the cylinder (0 < r < a), the expansion of the 
internal EM field (Ei,Hi) in terms of cylindrical har- 
monics M„ and N„ [lj| are expressed below. For the 
TM and TE modes, which are indicated by the indexes 
(I) and (II), respectively, one obtains 
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where E n = rE (— ^)' ^ / ( o 1 , with pi = kr(m 2 — cos 2 C)^ 2 > 
and the index (1) indicates the Bessel function J n {pi) 
to generate the cylindrical harmonics 19]. The quantity 
in = {piti/ pe) 1 / 2 is the relative refraction index between 
the cylinder and the surrounding medium, and k\ = mk 
is the wavenumber inside the cylinder. 

To simplify the expressions of the internal (c„, d n ) and 
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where £ = a;sin£, 77 = x(m 2 — cos 2 Q 1 ^ 2 , x = ka is the 

size parameter, Hn = Jn + iY n is the Hankel function, 
and fh — (/iei/ pit) 1 ! 2 is the relative impedance between 
the cylinder and the surrounding medium. These func- 
tions (fSI)— (fTUj) are analogue to those presented in [l9j], and 
they are the same for \i = pi (nonmagnetic approach). 

For the TM mode, the boundary conditions provides 
a set of four linear equations connecting the coefficients 
a ( n\ b ( n\ c ( n\ and d ( n ] [II [II S3- Solving the system 
of equations, we obtain: 
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Similarly, for the TE mode, we have: 
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where the functions A n , B n , C n , X>„, V„ and W„ are 
defined in Eqs. ([5])-([TD]). 

For normal EM wave incidence to the cylinder axis 
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(C = 90°), we have <$ = = ci n) = d^ = and [H 
_ mJ' n (x)J n (mx) — J n {x)J' n (m,x) 
C=9 °° ~ mJ n (mx)H^\x) - J' n (mx)H^ (xj 
J n (mx)J n '(x) -mJ n '(mx)J n (x) 
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J n (mx)H^ 1 \x) - mJ n ' (mx)H^\x) 
2i/ttx 

J n {mx)Hn 1 \x) — fnJ! n {mx)Hji\x) 
lijixx 



mJ n (mx)Hn l> (x) — J/ n (mx)H^ x> (x) 
where we have used the Wronskian Hn 1 \x)J n (x) — 

H [ n ) {x)J! n {x) = 2l/'KX. 

III. TIME-AVERAGED INTERNAL ENERGY 

The time-averaged EM energy within a nondispersive 
finite cylinder with radius a and length L is given by 
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This expression takes e± and \x\ as complex quantities 
with positive real parts and small imaginary parts com- 
pared to the real ones. In particular, for a cylinder with 
the same optical properties as the surrounding medium, 
one has: 

2 

W Q = ^-e\E \ 2 L. (20! 



To simplify the analytical expressions and, thereby, the 
numerical calculations, it is common to use some rela- 
tions involving the Bessel functions. Specially, for the 
average EM energy [25] provides two equations in which 
the integrals associated with product of two cylindri- 
cal Bessel functions are performed analytically. In our 
notation, for the situation in which there is absorption 
(m =^ to*), we can define the function 
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where p\(r) — kr(m 2 — cos 2 C) 1 ^ 2 an d f] = pi(a). Using 
the L'Hospital's rule and the recurrence relation J' n (p) = 
±[nJ n (p)/p — J n ±i(p)]j for real relative refractive index 
(to), Eq. ()21|) can be rewritten as 
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In addition, from the recurrence relations 2nJ n (p) = 
p[J n -l(p) +J n+ t(p)] and 2J' n (p) = J n -i(p) - J n +i(p), 
one can readily show that 
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for any functions A and B. Eqs. (J23J) and ([24|) are original 
and they appear in the calculation of the average energy 
associated with the components (r, <fi) of the EM field at 
oblique incidence. 

Consider the internal fields defined by Eqs. ([TJ-(|4]) and 
take separately each one of the field components in the 
definition (|T9")) . For the TM polarization, the average EM 

energy (a) is given by: 
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where W Er = / drr / d<j) f dzRe(ei)|£i r | 2 /4, W Hr = 
J drr J d(j) J dzRe(/xi)|iJi r | 2 /4, and so on, and D n {p\) = 

J'niPl)/ Jn{Pl)- 

Because of the integrals in the radial component, ob- 
serve that Eqs. (J2U), (@7]), ^ and dSOJ) cannot be solved 
analytically. However, if one considers the contributions 
W Er4> = {W Er + W E<P ) and W Hr * = (W Hr + W H<t> ) 
to the internal energy, the expressions can be simpli- 
fied by means of Eqs. ([2"3")l and ([2~4")l . Explicitly, using 

Eq. ([24]), for A = c^cosC/m and B = d$ , it follows 
from Eqs. (gU) and (JITJ) that: 



When m and m are real quantities, by using L'Hospital's 
rule, Eq. (f34|) takes the simple form 
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For the TE polarization (case II), we obtain similar 
expressions: 
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Also, employing Eq. (|2"3"|) . for A = c„ and _B 
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Note that the integral in the last term in Eq. ([3"2")l cannot 
be performed analytically. However, for normal incidence 
(C = 90°: Ei||*), it is clear that = W$ z = 0, 

and the average EM energy is W^ t — + Wjj r <k- 

In this particular case, it can be shown that T n ^i(y) + 
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For normal incidence (£ = 90°: Ej _L z), one obtains 
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If m and m are real quantities, Eq. (|4*U|) becomes Eq. ([55)1 
replacing c„ with g?„. 

In all equations above, we have used the equalities 
Re (mm) = Re (ex) /e, which is associated with the elec- 
tric field, and Re(mm*) = |fci/aj/ii| 2 Re (/ii) /e, which 
appears in the magnetic one. For dispersive cylinders, 
the expressions for the internal energy must be modified 
according to the model used to write the functions e\(uj) 
and mi (w) [H S3. 

IV. ABSORPTION EFFICIENCY 



and the coefficients Cn and dn^ , given by Eqs. (fl7|) and 
([18]). we obtain 
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where V„ is defined in Eq. ©. 

If we consider in the last term of Eq. (|46| the approx- 
imation rj 2 w r;* 2 , which means m 2 m* 2 (low absorp- 
tion), it can readily be shown that 
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The efficiencies in the EM scattering by a non-optically 
active infinite cylinder are 
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where <3to t , Qsca and Q® s are the extinction (or total), 
scattering and absorption efficiencies for the TM polar- 
ization, repectively. Expressions for the TE mode are 
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obtained replacing al with bn L> and bn 1 with an 1 ' [19j . 
Using the boundary conditions for the TM mode 
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and the definitions for the magnetic internal coefficients 
and dn\ given by Eqs. (fT5|) and (fl"4"|). we obtain 
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where W„ is defined in Eq. (TTU|) . In the same manner, 
using the boundary conditions for the TE mode 

vMv)^ = im$HW(0bM , (47) 
m^J^ry)^ 11 ) = m£ [j B (0 - H™®^] , (48) 



where we have used the Wronskian H£ l \£)J n (£) — 

Hn (QJniQ = 2«/7T^. An analogous expression is ob- 
tained for the TE polarization by replacing the index (I) 
with (II) in Eq. (|50l . 

It is important to emphasize that in Eq. (|50p only the 
terms that vanish for the normal incidence are approxi- 
mated by using rj 2 w rf 2 . Therefore, when £ = 90°, it 
follows that 
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which are exact expressions for the parallel and perpen- 
dicular absorption efficiencies expanded in terms of the 
magnetic internal coefficients c n and d n . In the fol- 
lowing, we show an expected interrelation between the 
average EM energy within a cylinder and its optical- 
absorption efficiency, provided that Re(m) 3> Im(m) and 
Re(m) 3> Im(m) (low absorption limit). 



V. NUMERICAL CALCULATIONS 

Here we present some numerical results from the ex- 
act expressions of the time-averaged EM energy within 
a magnetic cylinder. All numerical calculations have 
been performed by programs written for the free soft- 
ware for scientific computation Scilab 5.1.1. As an up- 
per limit N for the truncated series , , we employ 
the expression N = max(n c , \m\x) + (101.0 + a;) 1 / 2 , with 
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n c = x + 4.05X 1 / 3 + 2 [21]. The modification added in 
N in which we take the value max(n c , |m|a;) instead of 
only n c is introduced to give more accurate sums even 
for large values of Hi/fi at small values of x (Figs. [TJ [2] 
and [3J. The exceptions are the Figs. |H [5] and [6l where 
we replace max(rt c , |m|a;) with n c . 




Figure 1: Normalized EM energy W^/Wo within a magnetic 
(fii/fi = 10, 100) and a nonmagnetic (fJ,i/fJ> = 1) cylinder 
with ei/e = (1.334 + 1.5 x 10~ 9 ) 2 . Only the TM polarization 
is shown, with £ = 60°. 
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Figure 2: Comparison between the normalized EM energy 
Wtot/Wo within a magnetic (/J,i/fJ, = 100) and a nonmag- 
netic (mi//x = 1) cylinder with ei/e = 1.4161. The parallel 
and perpendicular polarizations are indicated by (I) and (II), 
respectively. 
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Figure 3: Comparison between the normalized energy- 
transport velocity ve/co in a medium containing magnetic 
(fii/fj, = 100) and nonmagnetic (fii/fi = 1) cylinders with 
ei/e = 1.4161 and volume fraction / = 0.36. The parallel 
and perpendicular polarizations are indicated by (I) and (II), 
respectively. 
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Figure 4: Normalized EM energy W^/Wq within a magnetic 
(fii/fj, = 10, 100) and a nonmagnetic (/ii//i = 1) cylinder 
with ei/e = 1.4161, parallel polarization. 

Fig. Q] shows a comparison between the magnetic and 
nonmagnetic approach for the TM mode with £ = 60° 
(oblique incidence). The quantities are calculated in the 
interval < x < 2, with Sx = 10~ 3 . We use in this 
calculation the same value of ei/e of [10] for a magnetic 
sphere with a small imaginary part added. The result 
achieved in Fig. [T] is quite similar to the one we have 
obtained for EM scattering by a magnetic sphere 
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Figure 5: Normalized energy-transport velocity ve/cq in a 
medium containing magnetic {fii/fJ, = 10, 100) and nonmag- 
netic (fii/fJ, = 1) cylinders with ei/e = 1.4161 and volume 
fraction / = 0.36, parallel polarization. 
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Figure 6: Normalized energy-transport velocity ve/co in a 
medium containing magnetic cylinders (fii/fi = 1000), paral- 
lel polarization. 



In both cases, the average internal energy is much larger 
than that one related to a scatterer with the same optical 
properties as the surrounding medium even for x < 1. 

The series of sharp peaks in Fig. U as a function of 
the size parameter are well-known and are generally re- 
ferred to as morphology-dependent resonances (MDR) 
[26| . For the average internal energy, they are ascribed 
to the resonances of the far-field scattering coefficients 
[27l . |28| , which are related to the internal coefficients by 



Eqs. (|44j), (|45[). (|47|) and fl48J). Physically, these large val- 
ues of the normalized internal energy can be explained by 
the enhancement of the extinction efficiency Qtot , whose 
sharp peaks occur at the corresponding size parameters 
in which the large values of W to t/Wo take place [l7| . 



A. Weak absorption 

Similarly to the case studied in [l2|,[2(|, it can be shown 
that for a cylindrical scatterer with weak absorption (wa) 
there is a relation between the average EM energy and the 
optical-absorption efficiency. For sake of simplicity, con- 
sider Eq. and assume that to, <C m r and to, <C fh r , 
where m = mr + trrii, fh = fh r + trhi are the complex rela- 
tive refraction and impedance indexes, respectively. Ap- 
proximating (to 2 — m* 2 ) ?» 4zTO r m.i, Re(mTO*) « m r m r 



and Re(mTO) 
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Comparing Eqs. (fSTj) and (|53l) . it results that 
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The relation (|54|) for cylindrical scatterers, which to the 
best of our knowledge has not been determined so far, 
holds for both TM and TE modes and it does not depend 
on the incidence angle (90° — £) . 

An analogous result has been obtained for the EM scat- 
tering by a dielectric sphere [2(3] and a magnetic one [l2j . 
Due to the system symmetry, the constant which asso- 
ciates Wtot/Wo with Qabs for a homogeneous magnetic 
sphere is not the same as that one for an infinitely long 



cylinder: W^ /W 



3m r Q^^ /(Srriix). In the spe- 



cial case of the infinite cylinder with optical properties 
similar to the surrounding medium, i.e., m r w 1 and 



W, 



(wa) 



Wq, it follows that: Q&bs 



nmiX. 



This last 



result is in agreement with [22|. For the sphere, one has 

Qab p s h) « 8m lX /3 [a nam. 

As a curiosity, comparing the EM scattering by spheres 
with infinite cylinders, both in the weak absorption 
regime (nii <C m r ), we can rewrite Eq. (|54l) as 



I abs 



Wo 



2V 



x rrii 

0"abs 



2kV mi 



(55) 



where V is the volume of the scatterer, k is the wavenum- 
ber of the incident EM wave, and a g and er a bs = Qabsfg 
are the geometrical and the absorption cross sections, 
respectively. Note that, although we explicitly consider 
two particular geometries in this derivation (a sphere of 
radius a and a segment L of an infinite cylinder of radius 
a), Eq. (|5"5)) does not depend on the shape of the scat- 
terer. Of course, Eq. ([5"5j) must be further investigated 
to verify if its universality is valid or not. 
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B. Energy-transport velocity 

Van Tigglen et al. [l8| have shown that, for simple 
3D dielectric scatterers, the energy-transport velocity Ve 
is related to the energy-enhancement factor Wtat/Wo by 
the expression ve = c / [1 + f{Wt t/Wo — 1)], where Co 
is the wave velocity in the host medium (e,/x) and / is 
the volume fraction occupied by the scatterers. Rup- 
pin [T3| has successfully used this expression for ii£ as 
an application of the average energy stored inside an in- 
finite dielectric cylinder. As it has been reported in (l7j . 
this simple model to calculate the transport velocity in 
2D medium reproduces well the results of Busch et al. 
]'2'.j . obtained in a different context of the low-density 
approximation of the Bethe-Salpeter equation. 

Here, we extend the use of this expression for ve to 
the calculation of the energy-transport velocity in a 2D 
disordered magnetic medium (£ = 90°). We consider a 
random collection of parallel isotropic cylinders with a 
packing fraction / = 0.36, which is the same used in 
the experiments with nonmagnetic scatterers Ti02 fl8j . 
Specially, we assume the scatterers are magnetic and have 
negligible losses, i.e., (ex,/ii) are real quantities. This 
last assumption can be achieved in soft ferrites, which 
present large values of with low magnetic losses at 
microwave frequencies typically below 100 MHz 

The quantities in Figs. [2] and [3] are calculated in the in- 
terval < x < 2, with 6x = 10~ 3 . The resonance peaks 
in the EM internal energy (Fig. [2]), as expected, pro- 
vide small values of the energy-transport velocity (FigJ3]) 
[]~7L [l8j . Because of the magnetism, the transport veloc- 
ity vanish even for cylinders with radius much smaller 
than the wavelength (Rayleigh size region). Vanishing 
of energy-transport velocity for x < 1 in a 3D disordered 
magnetic medium has been reported in Q ■ Qualitatively, 
this means that the EM wave spends a long time (dwell 
time) inside the scatterers, leading to a decrease in ve 
and, thereby, in the diffusion coefficient D = yp£*/3, 
where I* is the transport mean free path 0, 0, [29| . 
This strong decrease in the transport velocity, and con- 
sequently in the diffusion coefficient, is related to the 
single scatterer resonances and leads to electromagnetic 
wave localization [TH, [29J. The decrease of ve in a 2D 
disordered magnetic medium can be observed in Figs. [3] 
(for < x < 2) and [5] (for < x < 100). 

Notice that the energy-transport velocity plotted as a 
function of the size parameter shows an oscillatory be- 
havior. Because of the differences in the average EM 
energy between the parallel and perpendicular polariza- 
tions (Fig. [2]), the behavior of the normalized transport 
velocity is also different for both polarizations. Indeed, 
they show opposite oscillatory tendencies for small size 
parameters: while the oscillation amplitude of ve/cq in 
the parallel mode is reduced with increasing x, in the 
perpendicular mode it is increased. This can be clearly 
observed in Fig. [31 

In the interval 1 < x < 100, with Sx = 0.05, we 
show in Fig. [5] the profile of ve/cq (parallel mode) for 



Hi/fj, = 1000. Although the number of shaper drops 
in this size parameter region, for 8x < 0.05, is much 
larger than it is represented here, one can observe the 
global oscillatory behavior of Ve as a functions of x. In 
a particular configuration for x = 0.01, [7J has studied 
the quantity Ve/cq as a function of the relative magnetic 
permeability \x. 



VI. CONCLUSION 

The time-averaged EM energy inside an irradiated 
magnetic cylinder, for a general case of oblique inci- 
dence, has analytically been calculated for the TM and 
TE modes. We have shown that, similar to the dielectric 
[20] | and magnetic 12 1 spheres with low absorption, the 



optical-absorption efficiency associated with a magnetic 
cylinder, taken to be weakly absorptive, is related to its 
internal energy-enhancement factor. Indeed, this partic- 
ular result, when compared to that one from single Mie 
scattering, suggests a more general relation which does 
not depend on the shape of the scattering center. If its 
universality is valid, it can be applied to calculate in a 
simple way the energy-transport velocity in a disordered 
weakly absorptive media, provided only the volume and 
the absorption cross section of the particles embedded 
in the medium. Finally, we have used the normalized 
internal EM energy to determine the energy-transport 
velocity for 2D disordered magnetic medium. We have 
shown, in particular, the vanish of the transport velocity 
even in the Rayleigh size region. 
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Appendix A: Approximations for the normal 
incidence 



Below we consider the scattering coefficients in the 
small magnetic particle limit. Two other limits obtained 
are the Rayleigh and ferromagnetic ones. The scattering 
amplitudes are then calculated to address polarization in 
the Rayleigh limit. 



1. Small-particle limit 

If a; <C 1, we can approximate the values of the scat- 
tering coefficients a n and b n using the limiting cases of 
the cylindrical Bessel functions [19| • For simplicity, con- 
sider the cylinder is normally illuminated (C = 90°). The 
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values of these coefficients for n — and n = 1 are: 
b 



£± — -x (4 - x 2 ) + —x 2 (8 - x 2 ) 

8m v ' 32 1 ; 

iirmDo(mx) 



(4 - x 2 ) + gx 2 (8 - x 2 ) , 



-%irx 2 [m (8 - 3x 2 ) - D 1 (mx)x (8 - x 2 )] ( A1 ) 

32 [D\(mx)x + to] 
-ittx 2 [8 — 3x 2 — fhDi(mx)x (8 — x 2 )] 
32 [m,Di(mx)x + f] 



At = 90°, the value of degree of polarization P of the 
scattered light is not identically 1 for a magnetic cylinder: 



P = 



\mm — 1| 2 — |m/m — 1| 2 
\mrh — 1| 2 + \m/fh — 1| 2 



(A3) 



Therefore, similarly to the EM scattering by a magnetic 
sphere Q , the radiant intensity scattered by a magnetic 
cylinder is not symmetrical about 90°. 



2. Rayleigh and ferromagnetic limits 



In the Rayleigh approximation, x <C 1 and \m\x <C 
1, which leads to Do{mx) « —mx/2 — m 3 x 3 /16 and 
Di(mx) w 1/mx — mx/4, we obtain: 



a 



bo 



4 \m 



- 1 



i7r / m m 
32 \ m m 



— (mm — 1) x ( f + m m — 2mm) x 

4 v ; 32 v ; 



ai « - — 



6i — 



mm — 1 

TOTO + 1 

m/rh — 1 
m/m + 1 



?.7T 

32 
~ 32 



2m — 3mm + 1 



2m 



mm 

2 



1 

3m/m 



m/m + 1 



Note that, in the magnetic approach, the terms of order 
x 2 in the coefficients op and b\ do not vanish. The ex- 
pressions presented in [19J for the nonmagnetic scattering 
are recovered when m — rh. 

The ferromagnetic limit is also derived from the small- 
particle limit. However, we must consider both x <C 1 
and m|x 3> 1. For large arguments, the logarithmic 
derivative function can be written as 



D n (mx) 



1 



tan [mx - (2n + 1)tt/4] , (A2) 



and the ferromagnetic approximation is obtained by 
substitution of Do(mx) and D\{mx) calculated from 
Eq. (|A2j) into the far-field scattering coefficients of small- 
particle limit, Eqs. (|A1[) . 



3. Polarization 



For a normally irradiated cylinder, the degree of po- 
larization of the scattered light is P = Tu/Ti2, where 
T n = (ITxI 2 - |T 2 | 2 )/2, T 12 = (iTil 2 + |T 2 | 2 )/2, and the 
amplitude functions are T\ = i*o+2£™ =1 b n cos(n€>) and 
T 2 = a + 2 £~ a„ cos(n6), with 8 = 180° - <j> pj|. 

Using Eqs. (lAf [) in the Rayleigh limit, one can write 
the amplitude functions T\ and T 2 in terms of order x 2 : 

Ti w — - mm - f)x 2 - — —t— x 2 cos 9 , 

4 v ; 2 Vm/m + 1 / 

mm — I 

. m / i! \ mm 



„ m/m \ n vk mm — I \ „ 

r 2 «--(^-iU 2 -- — — — x 2 cos e 

4 V m / 2 V mm + 1 / 
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